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During this period I attended two Dynamics Explorer team meetings. I made a pre- 
sentation at the second team meeting on “How electric field ‘structures’ impact smooth 
F-region model calculations.” This presentation focussed on the problem faced by large 
scale ionospheric modelers who have to use smooth electric fields to simulate the real elec- 
tric fields which are rarely smooth. This work led to a paper which was submitted to 
the Journal of Geophysical Research. Also during this period I was able to extend the 
calibration of the DE 1 SAI images by working closely with Dr. M. H. Rees (Fairbanks, 
Alaska). I supported the work of Dr. J. U. Kozyra (Ann Arbor, Michigan) who is studying 
SAR arc phenomena using the DE 1 and DE 2 data base. 

Specifically, four objectives were accomplished during this reporting period. Three of 
these objectives were items associated with my DE proposal, while the fourth was a set of 
theoretical F-region computations supporting another DE investigator. 

(1) In the previous semi annual report it was stated that a first draft of a TDIM-DE 
study paper was distributed to co-authors. Work has continued along this avenue 
with meetings and correspondence between Dr. M. H. Rees in Fairbanks, Alaska 
and myself. We have compared the visible auroral image conversion algorithms 
used in this first study with those Dr. Rees has developed for data at different 
wavelengths. In the first study 630-nm and 557-nm images were used to deduce the 
auroral energy flux and characteristic energy of the precipitating auroral electrons. 
Dr. Rees’ new algorithms are applicable to various VUV filters and 557-nm. 
Indeed these new algorithms are needed for the second study period day 326, 
1981 since 630-nm is only available for two of the four DE 1 data sets. 

(2) The data for the second DE study (Southward IMF, B y negative conditions) was 
collected and put into global format. A total of 55 sets of auroral images were 
obtained from Dr. L. A. Frank’s group. Each set consisted of 3 images taken at 2 
visible and one VUV wavelength. These images were then converted to energy flux 
and characteristic energy data sets. Figures 1 and 2 show the auroral images at 
VUV (filter 123w) and 557-nm wavelengths, each is grey scaled and the intensity 
given in kR. This pair of images were obtained from the DE 1 SAI instrument on 
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the 21 November 1981 at 0354 UT. These data were combined using the algorithm 
of Dr. Rees to obtain the energy flux and characteristic energy, Figures 3 and 4, 
respectively. The previous method used to determine these two parameters was 
also used, since a 630 nm image was present. Figures 5 and 6 are respectively the 
energy flux and characteristic energy using this earlier method. There is general 
consistency between the two techniques, especially in the dynamic range of the 
energy flux (compare Figures 3 and 5). However the characteristic energies are 
somewhat different especially in the night sector auroral oval (compare Figures 4 
and 6). 

The convection pattern for this particular study period was put in a form compat- 
able with the global ionospheric model. A total of 32 plasma flux tube trajectories 
were selected to cover the high latitude ionosphere. Figure 7 shows these trajecto- 
ries in a polar MLT-magnetic latitude plot. This convection pattern was obtained 
from Dr. J. P. Heppner and corresponds to his southward, B y positive high lati- 
tude electric field model. His VEFI data set for D326, 1981 was found to be very 
consistent with this pattern. 

(3) The presentation made at the DE science team meeting on February 1987 was 
written up and submitted to the Journal of Geophysical Research. This paper 
discusses the shortcomings of representing the high latitude convection pattern 
as a smooth function. A structured electric field distribution was generated from 
published electric field observations. The consequence of such structures upon 
ionospheric densities was then studied. Attached is a copy of this paper. These 
electric field structures are being developed to help define the polar cap electric 
field distribution under B x northward conditions. 

(4) Dr. J. U. Kozyra requested theoretical information concerning the ionosphere’s 
dependance on SAR arc related phenomena. I carried out a series of mid latitude 
corotational model runs to quantitatively show how the F-region varied as a func- 
tion of electric field, topside number flux and a topside heat source. Attached is 
a copy of this brief report made available to Dr. Kozyra. 
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ABSTRACT 

A three-dimensional time-dependent ionospheric model was used to study how electric 
field structures affect the polar F-region. The electric field structures are represented 
by elongated Volland 2-cell models whose dimensions range from 10’s to 1000 km. These 
model structures are intended to represent the polar cap electric field for IMF B M northward 
conditions. A statistical method is used to generate a set of these structures. Their electric 
field strength and polarity are varied in order to study the F-region’s dependence on this 
magnetospheric input. For electric field structures whose size and electric field strengths are 
consistent with B a northward observations, several ionospheric dependencies were found. 
The N m F% parameter, in general, decreased due to the presence of the higher electric fields 
associated with the structures. Decreases ranging from a few percent to a factor of 4 were 
obtained from the simulations. This decrease in N m F% was further complicated by the 
initial N m Ft conditions, i.e., past history of the flux tube, being quite different when the 
structure is present. H m Fi in the vicinity of a structure can be raised or lowered by up to 
100 km, however, this change is not uniquely dependent upon the vertical induced drift. 
The role of enhanced reaction rates due to elevated ion temperatures is very important. 
The F-region becomes spatially very structured. The modulation of this structure in N m F 2 
is up to an order of magnitude in this study. 

1 . INTRODUCTION 

Our understanding of the individual physical processes at work in the ionosphere- 
thermosphere system is good. Both theoretically and experimentally the cause and ef- 
fect of many of these processes have been studied. Indeed, on a global scale numerical 
models allow the coupling of these processes to be studied, and the results compared with 
multi-parameter observations. However, on another level, the degree to which such models 
simulate the coupling is difficult to judge. Whenever the numerical model gives results not 
quite consistent with the observations, it is unclear as to whether the model inputs were 
i nsuffici ently constrained by the observations or whether the model’s spatial or temporal 
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resolution was adequate. Hence, in general, model-observation discrepancies can only be 
used with great difficulty and care to infer that the model's physical description is in- 
adequate. Therefore, determining how well the global ionospheric-thermospheric models 
reproduce the overall interactions is extremely difficult; a simple * error analysis ' is at best 
misleading because of its oversimplification. The objective of this study is to determine the 
global scale ionosphere’s response to the 'structure' in one of its inputs, i.e., the convection 
electric field. To date, no quantitative studies have been made to determine the effect of 
'structure' in global scale ionospheric models. 

Since the first high latitude ionospheric model study by Knudaen et a/., [1977] several 
large scale high latitude ionospheric studies have been undertaken [ Watkins 1978; Sojka et 
al., 1981a; Quegan et al. t 1982]. Although the model studies have become more sophisti- 
cated over the years, they all follow a similar methodology. They all use globally smooth 
input functions to represent the neutral atmosphere; neutral wind; and magnetospheric 
convection, precipitation, and heat flux. Due to the limitations associated with computer 
time and memory as well as observational coverage of these parameters, this methodology 
is quite reasonable. However, it is well known that such smooth inputs are the exception 
rather than the rule. The models simulate the global system in an averaged sense, and 
the small scale lengths, on scales of a few to thousands of kilometers, are smoothed away. 
To date, the effects that such structuring will have on the global simulations has only 
been hand-wavingly considered. Each structure is itself a self-consistent physical entity 
whose detailed characteristics are unknown. The accumulated effect of many structures is 
difficult to deduce because of the many different time constants involved. Since this study 
focusses attention on the magnetospheric convection electric field, a few examples from 
this area will be used to highlight the degree of structure associated with the high latitude 
electric field. 

In the pioneering work of Heppner [1972, 1973, 1977], tins very question was extensively 
addressed. Heppner characterized the smooth electric field signatures using a large data 
base of Ogo 6 polar region electric field measurements. He found several basic patterns 
and showed how they related to the magnetic activity and the interplanetary magnetic 
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field orientation. However, he also pointed out that on these patterns varying degrees 
of smaller scale electric field structures were superimposed. Furthermore, because of the 
abundance and strengths of the smaller scale structures, as many as 50% of the polar cap 
crossings could not be characterized unambiguously by a smooth pattern. Figure 11 from 
Heppntr [1977] dramatically demonstrates this problem. More recently, Frank et at. [1986] 
carried out a multi-parameter study of the disturbed polar cap; namely, during periods of 
Theta aurora. Using data from various instruments on the Dynamics Explorer satellites, it 
was again found that the electric field, as well as other parameters, were highly structured 
across the polar cap. Figures 5 through 12 of the Frank et al. paper show structure 
dominated patterns for each orbit. In addition to these papers, many other in-situ and 
ground-based observations have been made of the electric field structure. 

Table 1 lists pertinent features associated with the structures from a variety of these 
published observations. The amplitude refers to the departure from the average electric 
field value while the widths are based on a positive-to-negative peak separation. Both 
parameters are purely subjective, no detailed numerical analysis has been carried out. 
From Table 1 it can be seen that the structure range in size from 10’s of kilometers, 
(the lowest resolution obtainable from these published figures) to 1000 km; their average 
size being a few hundred kilometers. Also, the electric field variation from the average 
value ranges from a few mV/m to over a hundred mV/m. For these examples typical 
electric field variations of 10 to 30 mV /m were found. Unfortunately, these in-situ satellite 
measurements do not give a direct indication of the structure*s 2-dimensional form or 
its duration. To some extent ground-based incoherent radars viewing the polar cap may 
in the future be able to address this aspect. High resolution electric field observations 
from the Sondrestrom incoherent scatter radar [ Kelly and Vickrey , 1984; Koffman and 
Wickwar , 1984] and the EISCAT radar system [van Eykcn et aL, 1984] demonstrate that 
these radars have the necessary resolution to observe the spatial and temporal properties 
of local electric field structures. 

For these electric field structures to be included in a global scale ionospheric simulation, 
it is necessary to have a complete spatial and temporal description of the electric field. 
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Unfortunately, no such description is available. Theoretically, our understanding of the 
processes that lead to these structures is vague, although it is realized that the structures 
themselves are a complex ‘auroral- voltage-current’ system \Kan and Burke, 1985]. On the 
smaller scale size (sub-kilometers), treating the plasma structures (irregularities) as plasma 
turbulence has been reasonably successful [Kintner and Styler , 1985]. In attempting to 
describe plasma structures in a global F region model based on the transport formulation, 
these sub-kilometer irregularities [ Kintner and Seyler, 1985] are precluded and, therefore, 
only structure on scale sizes greater than 10km should be considered. The transport 
formulation does not allow for plasma waves, and, hence, conditions that lead to wave 
generation cannot be treated. Given that neither an adequate theoretical or observational 
description of the structures is available, a semi-empirical one must be created to carry 
out the study. 

To create a semi-empirical model of a polar cap structure, the mutual behavior of the 
precipitation, electric field, currents and conductivities must be known. At this time no such 
data base has been studied, although data bases such as the Dynamics Explorer’s could well 
prove -adequate. Instead, for simplicity, only one parameter (the electric field structure) 
will be modeled. As a result the study will not be directly applicable to the ionosphere; 
however, an indication as to how the model ionosphere quantitatively responds to multiple 
electric field structures will be obtained. Follow up studies for the other parameters and 
their interdependence would gradually complete the picture. 

In Section 2 the semi-empirical model of a single electric field structure is developed 
and specific structure polar cap electric field models are described. The results of using 
these models as an input to the USU Time-Dependent Ionospheric Model (TDIM) are 
presented in Section 3. Section 4 summarizes these results and discusses their implications 
with respect to present day ‘smooth’ input simulations. 
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2. Electric Field Structure 


2.1. Model 

In this section a very simple and analytic representation of the polar cap electric field 
structures is developed. This representation is only weakly based upon theoretical under- 
standing of the structures. The electric field structure is assumed to be associated with 
a discrete current system into and out of the polar ionosphere. Such polar cap discrete 
structures are often found to be sun-aligned. Other polar cap features such as “diffuse 
auroral patches* will not be represented by this numerical simulation. Further the electric 
field structure is assumed to be associated with finite length parallel sun-aligned sheets of 
current. The electric field pattern, to zeroth order, is then assumed to be a small elon- 
gated Volland 2-cell pattern [ Volland , 1973, 1975]. In order to quantify this electric field 
pattern the key parameters are chosen to generate a scale size and electric field amplitude 
consistent with observations. 

A Volland 2-cell electric field pattern consists of two regions; a central uniform elec- 
tric field region which in this case is elliptical in shape, and surrounding this a second 
region where the electric field diminishes to zero at some outer boundary (again elliptical 
in shape). The boundary between these two regions is the location of the hypothetical 
Birkeland current systems. If the 2-cell pattern is assu me d to be elongated in the sun- 
aligned direction, then the Birkeland currents are maximized at the dawn-dusk region of 
the boundaries between the two regions. Appendix 1 shows how the analytic model is 
constructed. 

Each electric field structure is represented by one such 2-cell pattern. The pattern is 
defined by a set of parameters: {xo,yo,»'o,/in»*’ 2 ,^o}. These parameters represent the 
following properties: 
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is the polar cap location at the center of a structure, 
is the central region semi-minor axis, 
is the ratio of the central regions semi-major axis to its 
semi-minor axis. 

is the outermost region semi-minor axis. At 
this boundary the electric field diminishes to zero, 
is the electric potential along the semi-minor axis. 


Given these 6 parameters the electric field structure can be generated. In appendix 1 
the units associated with these 6 parameters are defined. 

At this point polar cap observations of electric field structures have to be used to restrict 
the values of the 6 parameters in order to keep the analytic models in reasonable agreement 
with observations. In the references cited in the introduction it is noted that the structured 
electric fields are predominantly associated with IMF B z changing sign, zero, or positive 
(northward). Under B g (southward) conditions and almost irrespective of the activity 
level the polar cap becomes very quiet and unstructured. Table 1 lists the various electric 
field data papers used in determining realistic ranges for the electric field structures. For 
each paper the peak to peak distance, and the range of absolute electric field strengths 
are tabulated. Since these data were all obtained by in-situ observation from polar cap 
satellite crossings, they do not give direct information on either the polar cap distribution 
(x 0 ,yo)’s or the elongation of the “2-cell” model parameter (/k)’s. Clearly each parameter 
is observed to vary from one structure to the next. Therefore, these data must be used to 
generate in a statistical manner the probable distribution or variation for each parameter. 
The computer model then uses the probability function together with a random number 
generator to define the various parameters for the electric field structures. The structure’s 
spatial locations (xo»lfo) 876 distributed evenly over the polar cap while the remaining four 
parameters are assumed to have a value and a Gaussian distribution. A Gaussian 
random number generator is used for these latter four parameters. 
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Figure la shows the location and dimension of a randomly generated electric field struc- 
ture. The structure is located in the pre-noon sector and has a dawn-dusk central region 
of 550 km while the sun aligned dimension is extended to 920 km. For this structure a 
total potential drop of 12 kV was selected. Figure 16 shows contours of the electrostatic 
potential for this electric field model. Contours are drawn at 1 kV intervals with dashed 
lines representing negative potentials. This pattern is a basic Volland 2-cell with iono- 
spheric F-region plasma being transported antis unward in the central region. Also shown 
in Figure 16 is the adopted polar cap (x, y) coordinate system. The positive x axis has 
its origin at the magnetic pole and lies along the mid-night meridian. For both axes the 
units are degrees of magnetic co-latitude. Similarly, the y axis lies along the 0600 MLT 
magnetic meridian. 

2.2. Adopted Electric Field Models 

* t 

For this pilot study the objective is to determine quantitatively to what extent the polar 
cap electric field structures affect the F-region plasma characteristics. To date global iono- 
spheric simulations have adopted electric field models that do not contain the small scale 
structures described in the introduction. Hence, for this study 2 electric field structure 
models were constructed. They differ by the overall strength of the electric fields associ- 
ated with their individual structures. A third model will be briefly discussed. It is the 
reverse polarity of the stronger of the two models. In this light previous global simulations 
represent a fourth case for which the electric fields associated with the structures are zero. 
Such a case will be the baseline against which the new simulations will be compared. 

The information derived from published polar cap electric field structures (i.e, the con- 
tents of table 1) were used to define the model parameters distributions. These parameters 
are listed in Table 2 for the first two cases alluded to above. Both cases have identical spa- 
tial properties, their individual structures have a mean central width of 2° with a Gaussian 
half width of 1°; a mean ratio for the central semi major to semi minor axis of 1.7 with a 
Ga ussian half width of 0.3; and finally a mean outer boundary ratio to the central semi 
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minor axis of 1.5 with a Gaussian half width of 0.5. 

The Case A model has electric potential drops across the structures central regions 
defined by a mean value of 5 kV and a Gaussian half-width of 2 kV, while case B has a 
12 kV and 4 kV mean and Gaussian half-width, respectively. In order to complete the 
electric field structure model, a total of 45 individual structures were randomly distributed 
over a 12° radius polar cap area about the magnetic pole. For both Case A and B the 
same distribution of parameters for the 45 structures was adopted; they only differ in their 
electric field strengths. 

The spatial distribution of these structures and their respective central regions is shown 
in Figure 2 in the MLT-magnetic latitude polar plot. There is considerable overlap between 
the structures. This will generate rather irregular electric field distributions which will help 
simulate the observed electric field structure. Inside each central region in Figure 2 there is 
a number, this corresponds for each structure to the total central region potential drop for 
the case B model. Figure 3 shows contours of electric potential in the polar MLT-magnetic 
latitude frame for zero, case A and case B structure models. For these panels the structure 
electric potentials are added to the reference global scale 2-cell high latitude convection 
pattern. Top panel shows the Kp = 3 (60 kV) global potential pattern. Contours are 
drawn at 6 kV intervals with negative potentials being represented by dashed contours. 
The adopted large scale 2-cell pattern has a 5° antisunward offset and exhibits a slight 
throat pattern in the noon sector \Hcclis et al 1982]. Such convection patterns have 
been extensively used in large scale ionospheric and thermospheric modelling [Sojka et al ., 
1981a; Roble et al., 1982; Sojka and Sehunk, 1985a]. In Figure 3 the corotational electric 
field has not been added into the patterns. 

The middle panel of Figure 3 shows the case A electric field structure model. Although 
the contours are no longer as smooth as those in the top panel, they do, however, still show 
predominantly the same overall polar cap pattern. The electric field structure introduced 
in this panel has only caused distortions in the contours on the 6 kV scale (the contour 
level). In the lower panel which shows the case B electric potential the situation is quite 
different. Now the pattern does not resemble that of the top panel. There are many 
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contour ‘islands*. The potential pattern departs significantly from the large scale 2-cell 
pattern in the top panel. Inside the polar region the electric fields associated with the 
contours are no longer dawn-dusk or uniform. 

From these panels it is still difficult to visualize or even compare these electric fields 
with observations. Figure 4 shows the dawn-dusk electric field for the zero (top panel), 
case A (middle panel), and case B (bottom panel) electric field structure models. The 
top panel shows the 2-cell large scale electric field which is uniform across the polar cap 
with a field strength of ~ 12 mV /m. For case A and B models the polar cap is highly 
structured. In case A the electric field strength has a maximum of 25 mV/m departure 
from the zero case polar cap electric field. The electric field departure for case B reaches 
70 mV /m for one structure although more typically the departure is on the order of 40 
mV/m. These electric field strengths are quite consistent with the data given in Table 1. 
Furthermore the spatial scale of these structures ranges from about 1 to 5 degrees (110 to 
550 km at F region altitudes) along the dawn-dusk meridian. The spatial scale sizes are 
also consistent with the information in Table 1. Unfortunately at this time no detailed 
observational information is available on the spatial scales in a direction orthogonal to that 
shown in Figure 4. The model assumes these dimensions to be extended by about a factor 
of 1.7 over those in the dawn-dusk direction. 

One of the major functions of the electric field in the F-region is to drive plasma con- 
vection. The convection paths followed by the plasma are “(E* x B)/B 2 ” drift paths; 
where the electric field E/ is the vector sum of the electric field described above with the 
corotational electric field, and B is the magnetic field. Figure 5 shows 8 selected plasma 
convection trajectories for the zero (top panel), case A (middle panel), and case B (bottom 
panel) electric field structure models summed with the corotational electric field. The 8 
trajectories have been selected such that they all end equally spaced along the ‘x = 16’ line 
at 80,000 seconds UT. This is most clearly seen in the top panel where all 8 trajectories 
end in a line. Each trajectory end point is 2 degrees apart along the x = 16 line of the 
coordinate system defined in Figure 16. Such a choice of end points is somewhat arbitrary, 
but it would correspond to points along an approximate satellite track or lie inside the 
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field of view of an aximuth&lly scanning incoherent scatter radar. In both cases the data 
are all collected over a relatively short period of time, a few minutes at most. Having spa- 
tially and temporarily synchronised the trajectory end points, all earlier trajectory points 
are no longer synchronised between the three models. This being a direct consequence of 
the three different electric field models leading to different convection paths and speeds. 
For each of the three cases the trajectories are started outside the polar cap and where 
possible in sunlight. Hie final UT of 80,000 seconds was chosen because under the Kp = 3, 
2 cell convection pattern, top panel Figure 5, the starting points are all in sunlight giving 
reasonably uniform starting density conditions. Once these trajectories enter the polar 
cap they also enter darkness, the conditions which are most sensitive to stimulated density 
changes caused by the "(E* x B )/B 2 " drift mechanisms. Tick marks along the trajectory 
indicate hourly intervals. The trajectories take ~ 2 j hours to convect across the polar cap. 

In the middle and lower panels of Figure 5, the net polar cap convection is still anti- 
sunward although there are marked regions of sunward convection. The tick marks drawn 
at hourly UT intervals indicate that the transpolar cap times are changing significantly. 
Indeed this change has caused some of the day-side start locations to lie in darkness (due 
to the change in UT in this region). To avoid this the trajectory has been followed to 
earlier times when it was located in the morning sector auroral oval. Trajectory 8 suffers 
from this problem in both the case A and B model. 

The transpolar convection time differences are highlighted in Table 3, in which the time 
each trajectory takes to go from the x = —5.5 location to the x = 15.5 (« end point) is 
tabulated. These two locations are ±10° from the center of the unstructured convection 
pattern, which has an anti-sunward offset of +5°. For the zero structured case, the times 
range systematically from 1:38:18 hr:mn:s to 2:26:52. These increasing times represent the 
effect of the corotatkmal field in conjunction with the constant 2-cell polar cap electric 
field. Case A t imes vary randomly compared to the zero case, and with the exception 
of trajectory 1 and 7 are all shorter than the corresponding zero case trajectory. In the 
B case, the times are ail shorter, with trajectories 5 and 8 being about one quarter of 
the corresponding zero case times. The 30 min transport times correspond to an average 
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convection speed of 1.2 km/s. In going from no structures to case A the time changes 
appear larger than going from case A to case B even though for the latter comparison the 
structure electric fields have approximately doubled. This is a consequence of both the case 
A and B structures having the same area, and hence, the overall effect of changing just the 
potential is not as significant as changing the area and the potential. The modified transit 
times not only depend on the number and magnitude of the electric field structures, but 
also upon the area occupied by the structures. These faster transport times are expected 
to have significant implications for the polar cap density distributions. Naively, one could 
be tempted to infer higher densities along the faster moving trajectories because inside the 
polar cap the plasma is in darkness and decaying (i.e., reducing the time spent in darkness 
reduces the amount of decay). However, as the next section will show, the situation is 
considerably more complex even in the absence of plasma production sources. 

3. Modeled Ionospheric Parameters 

The three electric fields' described in the previous section were the variable inputs for 
three runs of the USU Time-Dependent Ionospheric Model (TDIM). This model and its 
history are described in section 3.1. Additional input parameters must be specified before 
the TDIM can be run; these are defined in section 3.2. The results of running the TDIM 
for each of the three electric field cases are given in subsections 3.3. and 3.4. 

3.1. Ionosphere Model 

The ionospheric model was initially developed as a mid-latitude, multi-ion (NO + , Oj”, 
Nj, and 0 + ) model by Sehunk and Walker [1973]. The time-dependent ion continuity 
and momentum equations were solved as a function of altitude for a corotating plasma 
flux tube including diurnal variations and all relevant E and F region processes. This 
model was extended to include high latitude effects due to convection electric fields and 
particle precipitation fay Sehunk et al. [1975, 1976]. A amplified ion energy equation was 
»1«r> added, which was based on the assumption that local heating and cooling processes 
dominate (valid below 500 km). Flux tubes of plasma were followed as they moved in 


13 

response to convection electric fields. A further extension of the model to include the 
minor ions N + and He*, an updated photochemical scheme, and the MSIS atmospheric 
model is described in Sehunk and Raitt [1980]. 

The addition of plasma convection and particle precipitation models is described in 
Sojka et aL [1981a, 6]. More recently, the ionospheric model has been extended by Sehunk 
and Sojka [1982a] to include ion thermal conduction and diffusion-thermal heat flow, so 
that the ion temperature is now rigorously calculated at all altitudes between 120-1000 
km. The adopted ion energy equation and conductivities are those given by Conrad and 
Sehunk [1979]. 

3.2. Additional TDIM Inputs 

Even though for this study attention is focused on the convection models, a number 
of other parameters, which are inputs to the ionospheric model, need to be described. In 
solving the equations for the ions listed above, a lower (120 km) boundary condition of 
chemical equilibrium and an upper (800 km) boundary condition of zero flux were imposed. 
The solar EUV and neutral atmosphere (MSIS) were specified by choosing the following 
solar and seasonal conditions; 1981 day 321, solar maximum, winter. Specifically, an A p 
of 30 and F10.7 of 170 were adopted for the studies. In addition, a simplistic neutral wind 
flowing antisunward over the polar cap with a speed of 200 m/s was used. This wind 
pattern is described by Sojka et al., [1981 a, 6]. 

» 

At high latitudes the other main magnetospheric input is the auroral oval, which acts 
both as a plasma production and heat source. The adopted auroral pattern for the case 
of 2 cell convection (Kp « 3) was the conventional Spiro et al. [1982] auroral energy flux 
model for an AE value of 400. Strongest auroral precipitation is present in the night sector 
oval with energy fluxes reaching 1.6 erg cm"* s _I . These regions of high energy fluxes all 
lie outside of the polar cap region and, therefore, do not significantly affect this study. 

Using the ionospheric model, the various inputs, and the convection cases described in 
section 2.2, three specific ionospheric density studies were made. For each run the eight 
trajectories shown in Figure 5 were followed over the polar cap as described. Each study 
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kept the solar thermospheric and magnetospheric inputs fixed for the entire simulation 
period, except for solar and thermospheric UT variations. In the following sections the 
three studies will be contrasted for a variety of ionospheric parameters. 

3.2.1. N m Ft 

Electron density profiles were obtained from the TDIM every few seconds (typically 
5 to 30 seconds) along each trajectory. As stated earlier each trajectory begins where 
possible with the appropriate sunlight steady state profile. Displaying data from the eight 
trajectories for the three models is somewhat difficult since the spatial coordinates and 
UT are no longer readily correlated between the runs. In fact for the selected trajectories 
only the end points have a correlation, i.e., by choice they all end at 80,000 seconds UT 
at fixed locations. Because the aim of the study is to look quantitatively at the effects 
of the structured E fields on the ionospheric plasma parameters, the normal method of 
presenting the global data, i.e., 2-D polar contour or color plots of specific parameters at 
a fixed altitude, is not useful because of its qualitative nature. Instead the parameters will 
be plotted directly as a function of the trajectory x location across the polar cap. In doing 
so the details of how each structure in the electric field affects the ionosphere can be seen. 

Figure 6 shows how the F-region peak electron density (7V m F 2 ) varies along each tra- 
jectory for the zero (top panel), case A (middle panel), and case B (bottom panel) model 
simulations. The densities are plotted against the trajectory x location. This being the 
location along the noon-midnight meridian direction, see Figure 16 for the definition of 
the axis. In the top panel where the plasma convection is primarily antisunward and uni- 
form, (see Figure 5 top panel) the trajectories move from the most sunward location, —10, 
over the pole to a final night side location of +15.5. During this trans-polar crossing the 
plasma moves into darkness with no particle precipitation. Therefore the plasma density 
is expected to decrease due to dominance of plasma recombination processes. Indeed this 
is observed; over this distance N m F% is found to decrease by about a factor of two. The 
trajectories begin with a range of N m Fi values which correspond to the different solar 
conditions at which each trajectory begins (see Figure 5 top panel). To zeroth order this 
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spread is maintained as the plasma flux tubes convect over the polar cap. At a higher 
resolution it is found that the final relative densities change by a few 10’s of percent. 
This latter difference is attributed to the slightly different transpolar times for these eight 
trajectories (see Table 3). 

Once the electric field becomes structured the systematic trends discussed above disap- 
pear; middle and bottom panels show this graphically. In the middle panel which corre- 
sponds to N m Ft along the eight trajectories shown in Figure 5 (middle panel) the densities 
show more than a factor of 4 difference between each other. The relative difference from 
x — —10 to x = +15.5 does not remain constant, and the initial densities are not the same 
as in the top panel Although the plasma still converts in a net antisunward direction, 
there are places where some of the trajectories skirt round electric field structures, see 
Figure 5, middle panel. In Figure 6, middle panel, these loops which include sunward 
excursions are represented by the N m F j curve becoming multi valued as a function of z. 
The shaded area labeled A highlights such a region. 

Clearly, the differences in density at the final locations for Case A are significantly 
larger than those for the zero case. In all 8 trajectories the case A final densities are lower 
than or equal to the corresponding zero case densities. What are the sources of these 
differences? Two mechanisms are primarily responsible and both stem from the presence 
of the electric field structures. First, due to the modified transpolar convection speeds the 
trajectories all begin on the sunward side of the polar cap at different UT’s and spatial 
locations. This results in quite different initial densities. Indeed two of the trajectories do 
not reach —10 on the z axis, see Figure 5 middle panel. Furthermore, the trajectory labeled 
8, which corresponds to trajectory 8 no longer reaches sunlight on the dayside hence its 
initial density is completely different from trajectory 8 in the zero case (top panel). In fact, 
the initial density for this trajectory is based upon the assumption that the trajectory has 
just left the auroral oval and that a steady state at the edge of the oval is equivalent to 
its earlier history. Clearly, no verification of this is possible without an extensive auroral 
precipitation and convection analysis. The consequence is that such trajectories can have 
a wide range of initial conditions that would depend upon their past history in the dark 
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auroral oval. 

The second cause of the enhanced density decrease is more directly associated with the 
electric field structure. In skirting round an electric field structure the N m Fi value shows 
a marked decay; see shaded area in Figure 6 middle panel. Even more marked density 
reductions are to be seen in the lower panel, the B case where the electric field structures 
are stronger. In subsection 3.2.3. the details of this enhanced interaction will be studied. 

Since only 8 trajectories have been studied for each case, it is exceedingly difficult to 
deduce statistically significant trends. Because of the large number of parameters that are 
involved in defining the structure electric fields, and then using random number generators 
to individually select these parameters, this study is not extensive enough to be treated as 
a Monte Carlo simulation. Hence, this study can only indicate the type of trends associated 
with the structured electric fields. 

Table 4 highlights the differences in N m Ft at the trajectory end points for the 3 cases. 
From Figure 6 it is evident that N m Ft depends sensitively on the structured electric fields. 
For all 8 trajectories of case A, the densities at the end of the trajectories were lower than 
or equal to the corresponding zero case. With one exception (trajectory 1), the case B 
densities were lower than the corresponding case A densities. 

The flux tube transit times for these trajectories also show similar trends. Six of the 
case A trajectories are shorter than the corresponding zero case, while 7 of the case B 
trajectories are shorter than the corresponding case A trajectories. As an indication of 
the stochastic nature of the undersampled simulation, the exceptional trajectories in these 
two comparisons are not correlated, i.e., the slower case B trajectory is not the higher 
density trajectory. This is further highlighted by the final density differences in Table 4. 
The differences in final N m Ft between case A and the zero case range from being the same 
(trajectory 6), to bring 3.2 times lower. Likewise, the difference between case B and case 
A extends from bring 1.14 times larger to a factor of 4.3 times lower. However, when the 
B case is compared to the zeroth case for the 8 trajectories all are between a factor of 1.6 
to 6.5 times lower. For this limited study, these specific factors can only be used to give 
a zeroth-order indication of the N m F 2 trends associated with the multiple electric field 
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structures. 

3.2.2. H m Fi 

Associated with the ionospheric N m Ft is the peak altitude, hmFi. Figure 7 shows the 
variation of hmFt *s * function of x for the sero (top panel), case A (middle panel), and 
case B (bottom panel) electric field models. Each panel shows all 8 trajectories. The 
layout of this figure is the same as for Figure 6. All eight zero case trajectories have to 
within a few kilometers the same /i m fj variation, decreasing from ~ 360 to 310 km. For 
the case A (middle panel) and even more so for the case B (bottom panel) hmFt varies 
significantly between the 8 trajectories. In case B hmFi approaches 500 km for several of 
the trajectories. The hmFi spread for case A is about 50 km while for case B it extends to 
150 km. These large h m F 2 departures are associated with the more intense electric field 
structures in which a vertical induced drift (Vj|) leads to h m Fi being raised or lowered. 

3.2.3. Influence of £ 

Figure 8 shows N m Ft, h m F 2 , Vjj, and E x along trajectory 7 for the case A electric field 
study. These four parameters are plotted against time, in hours, from the start of this 
trajectory. In this manner the time spent by the plasma flux tube in each structure is 
evident. Two events have been highlighted by shading and are labeled A and B. These 
two events lead to major changes in h m Fi and to a lesser extent N m Fi. 

The two lower panels in Figure 8 show details of the electric field structures. E ± is 
a measure of the horizontal electric field strength which not only leads to the plasma 
convection speed but also ion heating via enhanced ion neutral collisions as the electric 
field increases [see Schunk et al. t 1975, 1976). The electric field peaks at ~ 100 mV/m 
in event A. For this trajectory the average electric field is about 35 mV/m, which is 
significantly larger than the zero case value of ~ 12 mV/m In the panel above the bottom 
panel, the vertical induced drift due to the electric field is shown, upward drifts being 
positive. These vertical drifts range from 80 m/s (upward) to —75 m/s (downward). The 
effect of the strong upward drifts in events A and B is clearly evident on h m Fi (third panel 
from the bottom). H m Fi rapidly increases as Vj| increases. Event A leads to an hmFi close 
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to 500 km. Following the rapid increase in both Vj| and /i m Fj the recovery is significantly 
slower for h m Fi\ 10’s of minutes compared to the rise time of a few minutes for V\\. 

N m Fz (top panel) does not show as marked a signature during events A and B. Initially 
in event A, N m Fz drops by about 25% but then recovers somewhat, while in event B N m Fi 
shows a negligible response. In contrast N m Ft shows the most marked decrease in the 20 
minute period prior to event A. During this period i\T m Fj is decreased by a factor of more 
than 2, also in this period the average electric field is ~ 50 mV/m. This larger electric 
field acting for ~ 20 minutes is able to maintain a high ion temperature which maintains a 
rapid conversion of 0 + — ► NO + and, consequently, a fast recombination rate (see Sehunk 
et al., 1975, 1976]. As a result of this enhanced chemical reaction the F-region topside 
profile is being modified by a process other than diffusion. Indeed in Figure 8 immediately 
prior to event A, h m Fi is increasing even although Vg is downwards. The compositional 
changes in this location are dominant over the induced vertical downward drift. 

Another very important point to be drawn from Figure 8 is how complex and apparently 
uncorrelated are the changes in N m F 2 and h m F 2 . This is the case for a very restricted 
theoretical simulation in which only one input is allowed to vary. The overall changes in 
N m F 2 do not reflect the high degree of structure in the electric field. N m F 2 and to a lesser 
extent h m F 2 appear to filter out the marked variations on the seconds to a few minutes 
time scale. If N m Fi is looked at with higher resolution, it does respond at the few percent 
level to some of these changes. N m F 2 can respond relatively rapidly to only a small degree 
via the energy equation, where the electron and ion temperature time constants are a few 
seconds and 10’s of seconds, respectively. 

3.2.4. Ion Composition 

As alluded to in the previous sub-section, the F-region ion composition is radically 
altered in the electric field structure. This composition change arises from the enhanced 
rate of the 0 + — ► NO* reaction in the region of elevated ion temperature, itself caused 
by the enhanced electric field, hi Figure 8 the shaded region referred to as event A shows 
variations of N m Ft and hmFi which are not simply described by transport and diffusion 
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processes. This event is looked at with higher resolution as a series of altitude profiles of 
both 0 + and the molecular ions in Figure 9. Both panels show five sets of profiles, the 
molecular profiles are labeled 1 to 5. Corresponding pairs of 0 + and molecular profiles 
have the same line style. Ion density in cm -3 is plotted logarithmically to the base 10 
against altitude from 120 to 800 km. 

In panel A the five purs of curves represent a time sequence spanning 95 seconds 
starting with profile 1 and ending with profile 5. During this time the electric field is 
rapidly increasing from 29 mV/m to 90 mV/m (see figure 8 event A). Paralleling this 
electric field increase is a corresponding ion temperature (T,) increase from 1350 to 5000° 
K at 800 km. This electric field has generated a 2* hot spot [Schunk and Sojka, 19826]. As 
a result, the reaction rate between 0 + and NO + is enhanced causing 0 + to be converted 
to NO + . Hence the molecular ion profile which is primarily NO + at high altitudes is 
rapidly enhanced. This enhancement is highlighted by the large shaded arrow. NO + at 
high altitudes (> 350 km) is enhanced by one to three orders of magnitude. At 800 km 
the molecular ion density has increased from 3 to ~ 100 cm -3 . Also the height of the 
molecular peak has increased from 240 to 305 km. During this time 0 + is also being 
modified causing the bottomside to be rapidly eroded and effectively lifted by 50 km, and 
the small arrow indicates this trend. The 0 + peak is also being eroded away which causes 
it to appear to be lifting. This lifting of the peak and reduction in density is further 
enhanced by the increased T,-, which increases the scale height. H m Ft is also enhanced by 
the strong upward induced drift in event A (see Figure 8). As the 0 + peak is raised, from 
380 to ~ 480 km, its density is decreased by about 25%. 

In panel B (Figure 9) the subsequent “recovery” part of event A is shown. The electric 
field decreases from 90 mV /m to 35 mV /m in a 220 second period. Again the curve labeled 
1 represents the start of this period while curve 5 corresponds to the end of this period. 
During this period the T{ and Vj| are decreasing leading to a lowering of the peak. As T{ 
decreases the scale height reduces and this leads to a slight enhancement of the 0 + peak, 
in fact in N m Fi. The molecular ions no longer experiencing the favorable production from 
0 + , recombine in order to get back to what is primarily chemical equilibrium appropriate 
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to the Ti. Consequently, molecular ion density decreases and the peak is lowered. However, 
the time constant is somewhat slower for this process of recombination compared to the 
rapid 0 + — ► NO + production process. In the 3 1/2 minute period the molecular density 
above 350 km has only decreased by ~ 1 1/2 orders of magnitude. The decrease takes 
about four times longer than the original increase. 

An additional feature associated with the electric field structure is shown in Figure 9. 
At high altitudes (800 km) the molecular ion density is normally very low, typically less 
than 10 cm -5 (see curve 1 in Figure 9 panel A). However, during the event, in this case 
a period of about 4 minutes for the plasma flux tube to move through this structure, the 
molecular ion density is considerably enhanced. The electric field structure acts as a pump 
to drive molecular ions out of the bottomside ionosphere. It also acts as a heater. The 
topside ion temperature is raised to ~ 5000 °K in event A shown in Figures 8 and 9. 

3.3. Polarity of the Electric Field Structures 

t 

In setting up the electric field structure model, all 45 individual structures in both 
cases A and B were given a dawn-dusk electric field in their central regions. This choice 
produced antisunward convection through the center of each cell. Alternatively, it can 
be argued that during periods of northward IMF these polar cap structures give rise to 
sunward convection in their central region. Since no clear understanding of these structures 
is currently available, this latter situation would seem just as plausible as the assumed 
model of a structure. The question then is how sensitive is this study to the polarity of 
the electric field structures. To study this a further model case was computed. This was a 
repeat of case B with all 45 electric field structures having a dusk-dawn orientated central 
region electric field. For this reverse case, the structures central regions would tend to 
cause the plasma to convect sunward. The magnitudes of these electric fields would be 
identical to those shown in Figure 4, bottom panel. However, in this case the electric field 
polarity would be reversed, about a value of —12 mVm -1 , which is the electric field in the 
polar cap for the 2 cell pattern. 

Figure 10 shows the variation of N m Fi (top panel) and hmFt (bottom panel) for all 
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eight trajectories in the same format as the case B N m Fi and hmFt shown in Figures 6 
and 7 (bottom panels), respectively. Although each trajectory leads to a different value 
of N m Fi compared with the case B model, both have the same overall trend. The end 
point densities vary over the same range, reduced by a few percent to a factor ~ 8 of the 
density expected in the zero case. H m Ft undergoes changes which are also similar to those 
described for case B. It is enhanced up to almost 500 km. Overall, the result of reversing 
the polarity of the electric field structures is not significantly different to the plasma trends 
obtained for the original dawn-dusk orientated electric field structures. Specific details such 
as trajectory start locations, the plasma path through a structure, or where strong electric 
fields are experienced are quite different. 

4. Summary 

In this paper the F -region's response (as represented by the TDIM) to a series of electric 
field structures was studied theoretically. The individual electric field structures were 
modeled as elongated 2-cell Volland patterns. Their dimensions were small compared to 
the large scale high latitude convection pattern. Four sets of electric field structures were 
used to generate convection models: 

(a) Zero case, structures of zero electric field which corresponds to the conventional 
2 cell convection pattern. 

(b) Case A, a total of 45 structures whose mean electric potential drop was 5 kV. 

(c) Case B, the identical 45 structures in case A but whose mean electric potential 
drop was increased to 12 kV. 

(d) Reverse polarity case, identical to case B but with the polarity of the electric field 
in the central region being changed from dawn-dusk to dusk-dawn. 

For each of these four convection models the USU-TDIM was run for a group of 8 
transpolar plasma flux tubes. The plasma parameters associated with these runs were 
contrasted and the following results were obtained: 

(l) The overall convection process with the inclusion of the electric field structures 
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leads to different trajectories and transpolar crossing times. This implies that 
prior to entering the polar cap the initial plasma characteristics are different. 
Alone this difference can produce differences of factors of 4 in the overall N m Ft 
values. Initial conditions depend heavily on UT and location outside the polar cap 
especially during winter and equinox when the terminator’s location is important. 
Consequently, an understanding of plasma past history outside the polar cap 
would also be critical. 

(2) With enhanced electric fields being present in the polar cap for models case A, case 
B, and case B reverse polarity, the plasma polar cap transit times are generally 
reduced. Contrary to simple argument this faster transit time for a dark polar 
cap (no plasma production sources) does not yield higher densities. 

(3) The electric field structures lead to enhanced T,-, which in turn lead to enhanced 

0+ loss via the reaction 0 + — ► N0 + . This process is the dominant N m Fi 
mechanism. For the case B model it is responsible for factors of 4 lower s 

after the polar cap traversal. 

(4) The final N m Ft values for cases A and B are generally lower than the zeroth case. 
For case A the final densities seem randomly lower than the corresponding zero 
case. In contrast, for case B the final densities appear to be systematically lower 
than the corresponding zero case (see Table 4). This difference between case A 
and case B may not be important. Follow up studies will have to explore whether 
this difference is statistically meaningful. 

( 5 ) H m Fi experiences a variety of responses to both induced vertical drifts and rapid 
changes in composition and temperature. However h m Fi and N m Fz are not sys- 
tematically correlated, hi general, h m F^ is increased in the electric field structures 
by 10’s to 100’s of km depending on the electric field strength. 

(6) The JF-region peak and topside respond relatively slowly to small scale structures 
and, indeed, N m F% responds by only a few percent changes to sudden changes in 
the electric field. The integrated effect of the larger electric field (i.e., time spent 
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in such structures) leads to systematic decrease in the F-layer. 

(7) In contrast the topside molecular ion density, specifically NO + , responds dra- 
matically to sudden increases in the electric field. The NO* density at 800 km 
can increase by over 3 orders of magnitude in 10’s of seconds. Tins phenomena 
successfully pumps molecular ions into the bottom of the magnetosphere. 

Because of the large number of parameters used to define the structures, and then using 
random number generators to select these parameters, this study of 8 trajectories per case 
is not extensive enough to be treated as a Monte Carlo simulation. Hence, the results 
cannot be treated as statistically significant, but are merely indicative of the stochastic 
trends associated with the electric field structures. 

The overall impact of these electric field structures on the polar ionosphere is probably 
non-negligible. These structures exist during periods of IMF B„ = 0 or positive, however, 
how representative the model structures are, still needs further investigation. Without 
ascertaining tjie validity of the input structures, the conclusions can only give an indication 
of the possible types of ionospheric consequences. No theoretical description is available 
to extend the in-situ satellite data sets into a horizontal 2-dimensional picture or define 
their temporal characteristics. In this study, a two-dimensional two-cell model has been 
adopted with the assumption that these structures are stable for the length of time it 
takes a plasma flux tube to cross it. For the eight test trajectories, the average transit 
times through the individual structures range from 50 to 300 seconds. This study therefore 
assumes that the duration of an electric field structure is between 1 and 5 minutes. 

These test trajectories have been chosen such that they all follow different paths but 
arrive at the same time at their respective end points. If we had attempted to study 
the UT dependences of each trajectory, then the steady electric field structures would 
not be applicable. For studies in which a time evolution of the polar ionosphere is to be 
studied, the lifetime question of each structure must be addressed. Given that the structure 
model is not entirely unreasonable, it appears that the structures contribute to ionospheric 
processes in several ways; reduce N m f a, modify hmFt % locally enhance 2}, pump molecular 
ions to high altitudes, and lead to large-scale structuring in N m Fx. The degree to which 
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each of these ionosphere modifications depends upon the adopted structure model is very 
different. In describing the final N m F* variations, it is evident that two components affect 
it; first there is change in the initial conditions due to the different trajectories and transit 
times and, secondly, while in the vicinity of each electric field structure, drastic JV m Fj 
changes can occur. Hence, it is not only necessary to know the electric field details of 
each structure, but also to know how the structures are distributed and coupled. For 
example, by just changing the distribution geometry of the structures within the polar 
cap, which leads to different trajectories and, hence, transit times, would lead to modified 
initial densities and, consequently, final densities. Alternatively, keeping the distribution 
the same but changing the size or characteristics of the electric field would also lead to local 
changes in N m F * , causing the final densities to change. Other parameters such as T t or 
the composition, where the time constants are short compared to the trans-polar crossing 
times, would be more strongly dependent upon the details of the individual structures 
rather than the integrated effect of many structures. 

The above results are derived from a theoretical study of how electric field structures 
modify F-region plasma parameters. The study on its own is not a polar cap F-region 
study. It only indicates quantitatively the integrated effect of multiple electric field struc- 
tures upon the ionospheric N m F 2 , h m Fj and composition. These electric field structures 
are themselves only weakly representative of polar cap structures observed during periods 
of zero to northward IMF B x conditions. In order for this type of study to be directly 
related to the polar cap ionosphere, significantly improved understanding of the “structure 
mechanism” must be obtained. Schematically, Figure 11 shows the various processes which 
probably contribute to polar cap ionospheric structures. The strengths of the various cou- 
plings between the magnetospheric inputs are at this time relatively unknown. Several 
previous studies have been carried out to investigate how these other structure processes 
affect the .F-region; discrete auroral precipitation [Robinson et al. t 1985; Sojka and Schunk, 
1986]; neutral winds [Sojka and Schunk, 19855]. In the schematic, the shaded region in- 
dicates the sub-set of magnetospheric inputs which have been studied in tins paper for 
three adopted structure models. None of the couplings have been included. Significant 
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advancement in our understanding, either empirically or theoretically, would be required 
before such couplings can be included in this type of study. At the present time, extensive 
data bases, such as that of the Dynamics Explorer (DE) project, would be ideal for such 
empirical studies. Based upon the F-region modifications inferred from this study, it is 
very important that this type of modeling be extended to include these couplings with the 
other magnetospheric-ionospheric-thermospheric inputs. 
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APPENDIX 

The electric field structure Inside the polar cap Is defined as a sub-region whose center 
is at a MLT and magnetic latitude (AfLTo, DLATq). This location is taken to be the 
new origin. Cartesian coordinates are constructed as follows: +x direction is antisunward 
while the +y direction is dawnward parallel to the dawn meridian. A given location is then 
defined by (x p ,y p ) which can be transformed to a spherical polar latitude and longitude 
about this origin as follows: [x p and y p are in degrees of latitude at 300 km] 


lat = 90°-(*J + yJ)i 

(Al) 

long = tan~ 1 (y p /x p ) 

(A2) 


An Eulerian transformation can then be used to express this location in magnetic 
latitude-MLT coordinates. 

Centered on the displaced origin (xo>l/b) is an inner region of uniform electric field; this 
region being elliptical in shape. Outside tins inner region the electric field diminishes to 
zero at the outer boundary; again an ellipse. Inside the inner region the electric field is 
orientated along the dawn-dusk (or vice versa) direction. The inner region is extended in 
the noon-midnight direction. Figure la shows the orientation of these two boundaries. The 
inner region’s semi-minor axis (ro) lies parallel to the dawn-dusk direction. A parameter 
(/in) is defined to be the ratio of the semi-major (r m ) to semi-minor axis (f in = r m /ro). 
Inside this central region the electric field is constant and is defined by a total potential 
drop of <f > o kV. The electric potential 4>[x p ,y p ) is then chosen to be zero along the y p = 0 
line and hence; 

<H*p>Vp) = ( A3 ) 

An outer ellipse defines the region beyond which the electric field for this structure is 
zero. This outer ellipse has the same orientation and eccentricity as the inner ellipse. Its 


s emi- minor axis (rj) is given by rj = ro/ont- Any point lying between the two ellipses can 
be regarded as being at a radial distance r from the origin. This radius vector intersects 
the inner ellipse at rj and the outer ellipse at r n t< Equation A3 can then be used to 
determine the potential at and by definition the potential at r©* t = 0. Since the 
electric field diminishes linearly between rj and r***, the electric potential between these 
two boundaries 4>(r) is given by: 
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Table 1. In-situ observation of electric field structures. 


Reference 

Figure 

Satellite 

Amplitude (mV/m) Width (km) 

Heppner [1972] 

1 

Ogo 6 

10-20 

15-600 

Seppner [1977] 

11 

Ogo 6 

5-80 

10-800 

Burke et al. [1979] 

2 

S3-2 

10-20 

30-600 

HeeUt and Hanson [1980] 

4 

AE-C 

5-65 

few-1000 

Smiidy et al. [1980] 

5 

S3-2 

10-20 

40-400 

Frank et al. [1986] 

5-9 

DE-2 

10-100 

10-800 

Heelis et al. [1986] 

2 

DE-2 

15-60 

few hundred 
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Table 2. Adopted mean and half-width values for the 


electric field structures. 


Parameter 

Mean Value 

Half- Width 

central radius ro 

2° 

1° 

««»(a shs)* 

1.7 

0.3 

ratio rj/ro 

1.5 

0.5 

potential drop case A 

5 kV 

2 kV 

potential drop case B 

12 kV 

4 kV 
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Table 3. Convection times from x — —5.5 to * = 15.5. 


Trajectory 

Zero Case 
(Hr:mn:s) 

Case A 
(Hr:mn») 

Case B 
(Hnmna) 

1 

1:38:18 

1:49:00 

0:57:19 

2 

1:43:07 

1:33:03 

1:29:30 

3 

1:48:25 

1:33:23 

1:00:50 

4 

1:54:13 

1:22:20 

1:01:37 

5 

2:00:42 

1:03:35 

0:33:28 

6 

2:08:06 

1:40:55 

1:48:06 

7 

2:16:41 

2:42:22 

1:36:11 

8 

2:26:52 

1:01:00 

0:37:56 
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Table 4. Final N m Fj values*. 


Trajectory 

Zero Case 

Case A 

Case B 

1 

4.27 E5 

2.40 E5 

2.75 E5 

2 

5.37 E5 

3.89 £5 

3.16 E5 

3 

6.31 E5 

5.13 E5 

3.31 E5 

4 

5.89 E5 

5.25 E5 

2.69 E5 

5 

5.01 E5 

2.75 E5 

2.45 E5 

6 

4.68 E5 

4.68 E5 

1.07 E5 

7 

4.07 E5 

1.62 E5 

1.02 E5 


8 3.63 E5 

‘density given as logarithm base 10 (cm' 


f 
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Captions 

Fig. la. The area occupied by a modeled electric field structure (large ellipse) and its 
uniform central region (small ellipse) shown in a magnetic local time (MLT) and magnetic 
latitude frame. The number inside the small ellipse indicates the dawn-dusk potential 
across this ellipse in kV. 

Fig. 16. The corresponding electric potential contours drawn at 1 kV intervals shown in the 
new magnetic x-y coordinate system. Negative potentials are shown as dashed contours. 
The potential varies from — 6 to 6 kV for this electric field structure. The new coordinate 
system is located at the magnetic pole with the positive x axis lying along the magnetic 
midnight meridian while the positive y axis lies along the dawn magnetic meridian. Units 
along each aids correspond to degrees of magnetic colatitude. 

Fig. 2. Distribution of the adopted 45 polar cap electric field structures in the MLT- 
mngnetic latitude frame. Each structure is represented by its uniform central region and 
labeled in the center with its dawn-dusk potential drop in kV. 

Fig. 3. Polar cap potential contours of the unperturbed 2-cell convection pattern (top 
panel), case A model (middle panel), and case B model (lower panel) shown in the MLT- 
magnetic latitude frame. Contours are drawn at 6 kV intervals with negative potentials 
being represented by dashed contours. 

Fig. 4. Electric field variation along the dusk-dawn magnetic meridian for the unperturbed 
2-cell field (top panel), case A model (middle panel), and case B model (lower panel). The 
electric field is given in mV/m. 

Fig. 5. Plasma convection trajectories for 2-cell unperturbed study (top panel), for case 
A model (middle panel), and for case B model (lower panel) shown in the MLT-magnetic 
latitude frame. Circles denote the start location for each trajectory while tick marks 
indicate one hourly time intervals. 

Fig. 6. N m Fj densities for each of the trajectories shown in Figure 5 plotted as a function 
of the flux tubes’ x location. The panels correspond to the unperturbed 2-cell, case A 
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model, and case B model electric fields in the top, middle, and lower panels, respectively. 
The density is plotted as a logarithm to the base 10 of the density (cm “*). 

Fig. 7. H m Ft altitudes for each of the trajectories shown in Figure 5 plotted as a function 
of the flux tubes* x location. The panels correspond to the unperturbed 2-cell, case A 
model, and case B model electric fields in the top, middle, and lower panels, respectively. 

Fig. 8. Variation of the N m Fj (top panel), hm-Fs (upper middle panel), Induced upward 
drift (lower middle panel), and horizontal electric field magnitude (bottom panel) along 
trajectory 7 for the case A model electric field. The data are plotted as a function of 
elapsed time along this trajectory. Shading is used to highlight the two times when this 
trajectory passes through marked fluctuations in the induced vertical drift; the first event 
is labeled (A) and the second (B). 

Fig. 9. Variation of the 0 + and molecular ion density profiles for event A in Figure 8. 
The left panel corresponds to increasing disturbance while the right panel corresponds to 
the disturbance decaying. Unlabeled profiles represent 0 + while the profiles labeled 1 to 
5 are for the molecular ions. The shaded arrows indicate how the topside molecular ion 
abundance is changing with time, and the arrows shows how the bottom ride 0 + density 
varies with time. 

Fig. 10. N m Fi (top panel) and h m F^ (lower panel) variations along the 8 trajectories 
associated with the reverse case B model electric field plotted as a function of the flux 
tubes’ x position. The density is plotted as the logarithm base 10 of the density (cm -3 ). 

Fig. 11. Schematic of how structure mechanisms affect the F-region as represented by a 
physical model. The shading indicates the scope of this study. Dashed lines indicate links 
in the coupling processes which cannot be physically described at this time. 
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Center for Atmospheric and Space Sciences 
Telephone (801) 750-2961 
TELEX 3789426 


August 15, 1986 


Dr. Janet Kozyra 
Space Research Lab 
2455 Hayward 
Ann Arbor, MI 48109 

Dear Janet: 

It was nice meeting you again after so many years. I presume your return to 
Michigan was uneventful. The remainder of our holidays worked out well 
(getting rid of the kids for 3 weeks was great). 

Returning to USU regenerated and enthusiastic, I have followed up our COSPAR 
conversation concerning ionospheric signatures associated with SAR arcs. I 
have carried out a set of "trough" runs for different conditions, the 
differences being associated with your SAR arc signatures. I chose a total 
of eight trajectories, starting at noon local magnetic time, and followed 
them through midnight where they all arrive at 80,000 second UT. The runs 
are all for Solar Max (1981 F10.7 * 190), winter (Day 315) with Kp =* 5 and 
Ap « 40 which corresponds to active conditions. 

Figure 0 shows the trajectories. All the data plots correspond to 
latitudinal cuts just prior to midnight (the shaded area). Each of the runs 
is described separately but basically correspond to the following overview: 

Run Description 

1 normal high/mid latitude run with T 
frozen 

2 high/mid latitude run with T^ computed and 
top|ide^heat flux maximum of -3 x 10 eV 
cm s (downward flux) 

3 repeat of Run 2 with HF * -10*^ 

max 

4 repeat of Run 2 with (0]/lN-) ratio 
increased from 1 to 3 by [0] x 3.5 

5 repeat of Run 2 with convection = E x B + 
co- rotation + co- rotation 
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6 repeat of Run 2 with [N 1/3.5 to increase 
[Oj/lNj] ratio from 1 to 3 

7 repeat of Run 2 with convection ** E x B + 
co- rotation - 0.2 x co- rotation 

These runs roughly span the effects we talked about; they have relatively 
major effects on the topside ionosphere density* Please look over the plots 
and give me a call as to whether they make sense or not (801-750-2964) . 

The effects are significant enough that I would like to follow through with 
this study to produce a joint DE paper. As a DE guest investigator, this 
would be a rather appropriate study; how do you feel about it? Janet, since 
you made the suggestion and have the data, you will be the first author; at 
this end, I will be the only co-author. 


Best regards. 


Jan Sojka 
bw 

Enclosure 
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Figure 0. MLT - Magnetic Latitude plot of eight trajectories which are used 
for the studies. All begin at 12 MLT at different UTs but all 
enter the shaded region at 80,000 seconds UT. Study is Solar 
Maximum F10.7 * 190, 1981 with winter conditions (day 315); 
Kp = 5 and Ap *= 40. Tick marks are hourly intervals on each 
trajectory. Latitudinal plots which follow this are obtained in 
the shaded region. 
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Figure 2. Study 2 , a repeat <j>f the previous study but with T being solved 
and used in the 0 solutions. To solve at ^nigfitsid^, topside 
heat flux must be assumed a ^alue of -3 x 10 eV cm s was 

adopted. For this run the 0 topside (500 and 800 km) densities 
are increased by about 10%, not much. 
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Figure 3. Study 3 , same as ^tudy 2_^but^now the topside heat flux is 
increased to -1 x 10 eV cm s . 

The topside T has increased to 3000-4000° K. If the 
downward heat flux fs even "larger" then T will increase more. 

Mote, however, the 0 density increase at 500 km between 
Study 2 and 3 is^only about 10 to 30%. The largest T increases 
occur where 0 (hence n ) is smallest in the miSdle of the 
trough. 
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Figure 6. S tudy 6 , repeat of Study 2 and similar to Study 4, The [Oj/lN^J. 

ratio increased from 1 to 3 by tN^l/3.5. The topside 6 
densities are enhanced by about 30Z. 

Clearly, the importance of composition is evident, however, 
whether [0] is increased or [ 1 decreased would have to be 
determined empirically. 






DATA SET: 02 1 3K00M«SAR008.O 


RNRL41.F— (15 AUG 86) AT 12:02:43 


SEVENTH RUN (KP-5J . REPEAT RUN 2 BUT WITH CONVECTION- EXB • CORR -0.2«CORR 

300 KM 500 KM 800 KM 







MAGNETIC LATITUDE (DEG) 


Figure 7. Study 7 , repeat of Study 2 and similar to Study 5. The 
convection speed is changed by subtracting 0.2 x co-rotation 
speed. 

This small reduction in effective co-rotation speed results 
in a slight density decrease. The decrease being due to the 
longer time interval during which plasma decay occurs. 

For this and Study 5, no significant Joule heating effect on 
the ion temperature is evident. 





